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Abstract. This paper focuses on the study of open curves in a manifold M, and 
proposes a reparameterization invariant metric on the space of such paths. We 
use the square root velocity function (SRVF) introduced by Srivastava et al. in 
in to define a reparameterization invariant metric on the space of immer¬ 
sions M = lmm([0,1 ],M) by pullback of a metric on the tangent bundle T Jl 
derived from the Sasaki metric. We observe that such a natural choice of Rie- 
mannian metric on T M induces a first-order Sobolev metric on .M with an ex¬ 
tra term involving the origins, and leads to a distance which takes into account 
the distance between the origins and the distance between the SRV represen¬ 
tations of the curves. The geodesic equations for this metric are given, as well 
as an idea of how to compute the exponential map for observed trajectories in 
applications. This provides a generalized theoretical SRV framework for curves 
lying in a general manifold M. 


1 Introduction 

Computing distances between shapes of open or closed curves is of interest in many 
fields that require shape analysis, from medical imaging to video surveillance, to 
radar detection. While the shape of an organ or a human contour can be modeled by 
a closed plane curve, some applications require the manipulation of curves lying in 
a non flat manifold, such as S 2 -valued curves representing trajectories on the earth 
or curves in the space of hermitian positive definite matrices, where the values rep¬ 
resent covariance matrices of Gaussian processes. The shape space of planar curves 
has been widely studied (0,0, jT3) , 0), and the more general setting of shapes lying 
in any manifold M has recently met great interest ( (3j , [T2j, (5) , fl4l ) . Here we consider 
open oriented curves in a Riemannian manifold M, more precisely the space of im¬ 
mersions c: [0,1] —i• M, 

M - lmm([0,1 ],M). 

Reparameterizations will be represented by increasing diffeomorphisms <p : [0,1] —«• 
[0,1] (so that they preserve the end points of the curves), and their set is denoted 
by Diff l_ ([0,1]). Then, one way to describe a shape is as the equivalence class of all 
the curves that are identical modulo reparameterization, and the shape space as the 
associated quotient space, 


SP — lmm([0,1],M)/Difi^([0,1]). 


The formal principal bundle structure n : Jl —► SP induces a decomposition of the 
tangent bundle TJl — V Jl ® HJl into a vertical subspace V Jl = ker(T7r) consist¬ 
ing of all vectors tangent to the fibers of Jl over 5f, and a horizontal subspace HJl = 
(VJf) ±G defined as the orthogonal complement of V Jl according to the metric 
G that we put on Jl. We say formal because the manifold structure of the space 
lmm([0,1],M) has not yet been thoroughly studied to our knowledge. We require 
that G be reparameterization invariant, that is to say that the action of DifT 1 " ([0,1]) 
be isometric for G 

G C o<i>{h°(!),ko(p) = G c (h, k), (1) 

for any curve c £ Jl, reparameterization <p e Diff^QO, 1]), and infinitesimal defor¬ 
mations h, k e T C .M - h and k can also be seen as vector fields along the curve c in 
M. That way, the induced geodesic distance between two curves c 0 and c t does not 
change if we reparameterize them the same way, that is 

d{c 0 o(p, c\o<p) = d(co,Ci), 

for any <p e Diff + ([0,1]). Also, if this property is satisfied, then G induces a Rieman- 
nian metric G on the shape space, 

G 7 i(c) ( T c n{h ), T c ji{k)) = G c (h H , k H ), 

in the sense that the above expression does not depend on the choice of the repre¬ 
sentatives c, h and k. Here h H , k H denote the horizontal parts of h and k according 
to the previously mentioned decomposition, as well as the horizontal lifts of T c 7i{h) 
and T c n(k), respectively. The geodesic distances d on .M and d on J are then simply 
linked by 

d{[c n ], [cil) =inf{d(c o ,cio0) | 0eDiff^ftO, 1])}, 

where [co] and [ci] denote the shapes of two given curves cq and ci, and d verifies 
the stronger property 

d{CQ°<p, Ci o y/) = d(c 0 ,ci), 

for any reparameterizations (p,y/ e llifr([0,11). The most natural candidate for a 
reparameterization invariant metric G on Jl is the L 2 -metric with integration over 
arc length, but Michor and Mumford have shown in (6j that the induced metric G 
on the shape space always vanishes. This has motivated the study of Sobolev met¬ 
rics (ISMU.ED). and particularly of a first-order Sobolev metric on the space of plane 
curves, 

G c {h, k) - j (D e h- L ,Dfk- L ) + ^(D c h 11 .Dfk 11 ) (2) 

where we integrate according to arc length (if: = || c' ( t) || At and (■, ■) denotes the eu¬ 
clidean metric on IR 2 , D(h- h' is the derivation of h according to arc length, 
D s h 11 = < D s h , v)v is the projection of D s h on the unit length tangent vector field 
v ~ Wi^ a ^ on 8 c ’ anc ^ Dsh 1 - = < D s h , n) n is the projection of D s h on the unit length 
normal vector field n along c. This particular first-order Sobolev metric is of interest 
because it can be studied via the square root velocity (SRV) framework, introduced 
by Srivastava et al. in [lT| and used in several applications (0), (123). This framework 


can be extended to curves in a general manifold by using parallel transport, in a way 
which allows us to move the computations to the tangent plane to the origin of one 
of the two curves under comparison, see (5] and fl4l . In QD the transformation used 
is a generalization of the SRV function introduced by Bauer et al. in (T) as a tool to 
study a more general form of the Sobolev metric (2). In fT4] a Riemannian frame¬ 
work is given, including the associated Riemannian metric and the geodesic equa¬ 
tions. While our approach in this paper is similar, we feel that the distance we intro¬ 
duce here will be more directly dependent on the "relief" of the manifold, since it is 
computed in the manifold itself rather than in one tangent plane as in 0 and (T4]. 
This enables us to take into account a greater amount of information on the space 
separating two curves. 

2 New metric on the space of parameterized curves 

We consider the square root velocity function (SRVF) introduced in (TT) on the space 
of curves in M, 

c' 

R:M^TM, c >->■ _= , 

Vidi 

where || • || is the norm associated to the Riemannian metric on M. This function will 
allow us to define a metric G on M by pullback of a metric G on TM. First, we define 
the following projections from TTM to TM. Let £ e fp,u) TM and (x, U) be a curve 
in TM that passes through (/;, u) at time 0 at speed f Then we define the vertical and 
horizontal projections 

v P(p,u) : T(p,u)TM —• T p M, £ —>• £v:= V^o )U, 
hp (p,u) ' T(p,u)T M —*■ TpM, £x (0). 

The horizontal and vertical projections live in the tangent bundle TM and are not to 
be confused with the horizontal and vertical parts which live in the double tangent 
bundle TTM and will be denoted by % H , £, v . Furthermore, let us point out that the 
horizontal projection is simply the differential of the natural projection TM —► M, 
and that according to these definitions, the Sasaki metric (|9), fTOl l can be written 

8tp,u)G> T l) = (SH,TlH) + {Sv,Tlv), 

where (■, •> denotes the Riemannian metric on M. Now we can define the metric that 
we put on TM. Let us consider heTM and £,77 e T^TM. We define 

G h (dT]) = (CO)//, 77(0)/,) + f 1 ( 3 ) 

Jo 

where £(?)# = hp(£(f)) and £( t)y = vp(<f(f)) are the horizontal and vertical projec¬ 
tions of ( t) e TTM for all t. Then we have the following result. 

Proposition 1. The pullback of the metric G by the square root velocity function R is 
given by 


G c [h, k)-{h{0), k(0) > + J (\! c h L ,\J e k L ) + (h 11 ,V e k l,S ) AT 


(4) 



for any curve c e Jl and vectors h,lc e T c Jl, where we integrate according to arc 
length, Vf h = jpjV c 'h is the covariant derivative ofh according to arc length, and 
V (h^ = {D(h, v)v and -V(h-V(h^ are its tangential and normal components 
respectively. 

Remark 1. In the case of curves in a flat space, G is the first-order Sobolev metric (2), 
studied in ED , with an added term involving the origins. This extra term guaranties 
that the induced distance is always greater than the distance between the starting 
points of the curves in M. 

Proof. For any c e Jl, and h,ke T c Jl, the metric G is defined by 
G c (h, k) = Gr( C ) [T c R{h), T c R{k)). 


For any t e [0,1], we have T c R{h){t)H = h{i) and T c R[h)v = VhR(.c)[t). To prove this 
proposition, we just need to compute the latter. Let 5 >-► c[s, •) be a curve in Jl such 
that c(0, ■) - c and c s (0, ■) - h. Here and in all the paper we use the notations c s = 
dc/ds and c t — dcldt. Then 


Vfcfl(cXf) = 


1 


||c'|| 1/2 

1 

l|cdl 1/2 

1 

llc t || 1/2 


|| 1/2 


Vfcc' + h||| c'|| 1/2 )c' 

V s c t + d s (c t) c f > _1/4 c t 
V f c s -i(c t ,c t >- 5/4 (V s c t ,c t >Cf 


iV(h) ± + -{Vth, ->- 

2 W lie'll lie'll 


where in the last step we use again the inversion V s c f = V t c s . 


3 Fiber bundle structures 

Principal bundle over the shape space We already know that we have a formal prin¬ 
cipal bundle structure over the shape space 

7t : Jl — lmm([0,1],M) — & = JtlDift ([0,1]). 

l 

which induces a decomposition TJ{- V Jl ® HJl. fust as in the planar case, the 
fact that the square root velocity function R verifies the equivariance property 

R{c °</>) = [R{c) o (/)) 

for all c e Jl, h, k e T c Jl and <p e llifr ([0,1]), guaranties that the integral part of 
G is reparameterization invariant. Remembering that the reparameterizations <p e 
Diff 4 " ([0,1]) preserve the origins of the curves, we notice that G is constant along the 
fibers, as expressed in equation (l), and so there exists a Riemannian metric G on 
the shape space SP such that n is (formally) a Riemannian submersion from {,11., G) 
to (J, G) 

G c {h H , k H ) = G n(c) {T c ttfh), T c ttm, 
where h H and k H are the horizontal parts of h and k respectively. 





Fiber bundle over the starting points The special role that plays the starting point 
in the metric G induces another formal fiber bundle structure, where the base space 
is the manifold M, seen as the set of starting points of the curves, and the fibers are 
the set of curves with the same origin. The projection is then 

c^c(O). 

It induces another decomposition of the tangent bundle in vertical and horizontal 
bundles 


V^Jl = kerTTr'" 0 = {h£ T c M\h(Q) = 0}, 

H^Jl = (V C M ^) ±G . 


Proposition 2. We have the usual decomposition TM - V^M ® H^M, the hori¬ 
zontal bundle consists of parallel vector fields along c, and is (formally) a 

Riemannian submersion for (Ji, G) and ( M , (-, •>). 

Proof. Let h be a tangent vector. Consider ho the parallel vector field along c with ini¬ 
tial value ho(0) = /? (0). It is a horizontal vector, since its vanishing covariant deriva¬ 
tive along c assures that for any vertical vector l we have G c (/to, Z) = 0. The difference 
h-h - ho between those two horizontal vectors has initial value 0 and so it is a ver¬ 
tical vector, which gives a decomposition of h into a horizontal vector and a vertical 
vector. The definition of hf* ) ..M as the orthogonal complement of V { * ] .M guaranties 
that their sum is direct. Now if k is another tangent vector, then the scalar product 
between their horizontal parts is 

G c {h H , k H ) = ( h H ( 0), k H { 0) > c(0) = (MO), jfc(0)) c{0 ) = ( T C 71 M ( h H ), T c ti m ( k H ) ) nM , 
and this completes the proof. 


4 Induced distance on the space of curves 

ffere we will give an expression for the geodesic distance induced by the metric G. 
Let us consider two curves Cq, cy e M, and a path of curves s >-► c[s, •) linking them in 

M 

c(0, t) - co{t), c(l, t) = Ci(f), 

for all t e [0,1]. We denote by q{s, •) = R (c(s, •)) the image of this path of curves by the 
SRVF R. Note that q is a vector field along the surface c in M. Let now q be the raising 
of q in the tangent plane T C (o,o)M in the following way 

^ f ) = p c s ;°,o ) oP Z,-M s ’ r) )’ 


where we denote by Py ’ 12 : Ty^pM —>• 7j,( (2 ) M the parallel transport along a curve y 
from y(fi) to yltf). Notice that q is a surface in a vector space, as illustrated in Figure 



[I] Lastly, we introduce a vector field ( a , t) >-► of' 1 (a, r) in M, which parallel translates 
q{s, t) along c(s, ■) to its origin, then along c(-,0) and back down again, as shown in 
Figure[l] More precisely 


(D ' (a, T) “ * cf.,0) ” c(s,0 


,o P 


s,a 


o P 


t, 0 


[q(s, r)) 


for all b, s. That way the quantity V s cj s,t measures the holonomy along the rectangle 
of infinitesimal width shown in Figure^ We can now formulate our result. 

Proposition 3. With the above notations, the geodesic distance induced by the Rie- 
mannian metric G between two curves cq and c\ on the space = lmm(\(), 1 ],M) of 
parameterized curves is given by 


dist{co,c\) - inf f t/l|c s (s,0)|| 2 + f ||V s a(s, f)|| 2 dt ds, (5) 
c(0,-)=c 0 ,c(l,-)=ci Jo V Jo 

where q - R ( c) is the Square Root Velocity representation of the curve c and the norm 
is the one associated to the Riemannian metric on M. It can also be written 

dist(c Q ,C\)- inf f \ \\c s {s,0)\\ 2 + f II <5 s (s, t) + Q(s, t) II 2 df ds, (6) 

c(0,-)=c 0 ,c(l,-)=ci J 0 y J 0 11 11 

where q is the raising of q in the tangent plane r C (o,o)M and the curvature term Q is 
given by 

m*'t) = P s £ 0) oP c f ; 0 v) (V s ^ f (s, f )) 

= C ° P cL (f‘ P cL fa c *’ C ^ P c(l>^’ «) d ") - 

iff? denotes the curvature tensor of the manifold M. 


Remark 2. Our original motivation for this work was to find a geodesic distance (that 
is, a distance induced by a Riemannian metric) that resembled the product distance 
introduced in (5) . In the first term under the square root of expression (6) we can see 
the velocity vector of the curve c(-,0) linking the two origins, and in the second the 
velocity vector of the curve q linking the TSRVF-images of the curves - Transported 
Square Root Velocity Function, as introduced by Su et al. in fl2l . ffowever there is 
also a curvature term Q which, as previously mentionned, measures the holonomy 
along the rectangle of infinitesimal width shown in Figure[l] If instead we equip the 
tangent bundlewith the metric 


GfcGU) 




pyimc'.ShipyMt)) dr 


2 

dt, 


for h e TM and f i] e I), T.M , then the curvature term Q vanishes and the geodesic 
distance on Jl becomes 


dist(co.ci) 


inf 

c(0,-)=co,r(l,-)=ci 



llc s (s, 0) || 2 + ||<7 s (v)|| 2 ds, 


(7) 







where the norm of the second term under the square root is the L 2 -norm, and which 
corresponds exactly to the geodesic distance associated to the metric on the space 
C = U peM L 2 {[0, 1], TpM ) introduced by Zhang et al. in |[T4). The difference between 
the two distances {5} and {7} resides in the curvature term Q, which translates the fact 
that in the first one, we compute the distance in the manifold, whereas in the second, 
it is computed in the tangent space to one of the origins of the curves. Therefore, the 
first one takes more directly into account the "relief" of the manifold between the 
two curves under comparison. For example, if there is a "bump" between two curves 
in an otherwise relatively flat space, the second distance (7) might not see it, whereas 
the first one (5) will thanks to the curvature term. 

Remark 3. Let us briefly consider the flat case : if the manifold M is flat, the two 
distances (3 and {7) coincide. If two curves cq and c\ in a flat space have the same 
starting point p, the first summand under the square root vanishes and the distance 
becomes the L 2 -distance between the two SRV representations q () = W(c 0 ) and < 7 ! = 
R(ci). If two curves in a flat space differ only by a translation, then the distance is 
simply the distance between their origins. 


<w s '(s,0) 



c( 0,) = c 0 c(l,) = Ci 


Fig. 1. Illustration of the distance between two curves cq and c\ in the space of curves M 


Proof. Since G is defined by pullback of G by the SRVF R, we know that the lengths 
of c in M and of q - R(c) in 'X.M are equal and so that 


dist(c 0 ,ci) 


inf 

c(0,0=co,c(l,-)=ci 



(v),<7i(s.O) ds, 


G(<7s(vb<7s(v)) = ||c s (s,0)|| 2 + / ||V s ^(s, f)|| df. 


f I 


with 






Now let us fix t e [0,1], Then s >-► P*'® s f [q{s, f)) is a vector field along c(-,0), and so 




=p 


,0,s 

c(-,0) 


_ P S ’° o p 

c(-,0) c(,s 


C(5,0 S ’ “ P d; 0) i^ s ^ S ’ ’ 


We consider the vector field v along the surface (s, x) c(s, x) that is parallel along 
all curves c(s, ■) and takes value v(s, t ) = q(s, t) in x = f for any s e [0,1], that is 

v(s,x) = P*£ 0 (<7(s,0), 

for all se [0,1] and x e [0,1]. That way we know that 


V s v(s, t) = V s q{s,t), 
V s v(5,0) = P°’* 0) [q s {s, t )), 
V r v(s,x) = 0, 


for all s, x e [0,1]. Then we can express its covariant derivative in the following way 

5 ° ,f m .,.^J*p^^ tV ( s , T))dT 


V*vCs,r) = P“£ 0 (V s v(s,0)) + 


= P cL° P cm «) + [pl'l } (mc T ,c t )P%'. ) qts, t)) dx. 
«/ 0 


( 8 ) 


Now let us fix s e [0,1] as well. Notice that the vector field oj s,t defined above verifies 


cu s,t (s, t ) = q{s, t), 


V T w s,f (rt,x) = 0, 
V a w s ' f (fl, 0) = 0, 


for all a, x e [0,1], Note that unlike v, we do not have V s w s,t (s, t) = V s £/(5, t) because 
cj s,t [a, t ) = q{a, t) is only true for a = s. It is easy to verify that the last term of equation 
{§) is precisely the covariant derivative of the vector field to s,t 

V s w s,f (s, t) = f)] dx, 

since for any x e [0,1], a> s,t (s ,x) = P*’*t), and finally by composing by P^° 0) ° 
Pels .)> we obtain the second expression {6), which completes the proof. 


5 Geodesic equation on M 

In order to be able to compute the distance given by between two curves, we first 
need to compute the optimal deformation .v *-» c(s, •) from one to the other. In other 
words, we need to characterize the geodesics of .M for our metric. In order to do so, 
we use the variational principle. The beginning of the calculations are very similar 
to those in (Hi. Let us consider two curves cq, c\ e Jl and a path [0, l]3s« c(s, •) e 


.M. going from one to the other. This path c is a geodesic if it minimizes the energy 
functional E: M —«■ K + 



Let a !->• a e {-£,£), be a proper variation of the path 5 c(s, •), meaning that 

it coincides with c in a - 0, and it preserves its end points 


c(0, v ) = c, 

c a {a, 0 ,-) = 0 \/ae {-£,£), 
c a [a, 1,-) = 0 \/ae {-£,£). 

Then c is a geodesic of Ji if 4-\ E{c{a ,-,•)) = 0 for any proper variation a 

aa \a=0 

c{a, ■, •). If we denote by E(a) — E(c(a, •,•)), for a £ {-£,£), the energy of a proper vari¬ 
ation of c, then we have 


E{a) - J (c s (a,s,0),c s {a,s,0))ds + J j (V s q{s,t),V s q{s, t))dtds, 


where q = CtIVWctW is the SRVrepresentation of c. Its derivative is given by 

= J iV a c s (a,s,0),c s (a,s,0)) ds + f J {V„V s g(a,s, r),V s g(a,s, f)) drds. 
Considering that the variation preserves the end points, integration by parts gives 
/<V„c„c,>d s =-/<V,c„c„> ds 
J(V s V a q,V s q)ds = - J (V s V s q,V a q)ds, 


and so we obtain 


^E'(a) = - J <V s c s ,c fl )l r= (ids + J J (@(c a ,c s )q + V s V„d,V s d)dfds 

= -/<V,c,.c„>l,. 0 d S -//(« ( ,,V s „ C „c a > + (V ) V„,V a ,>d,d, 


This quantity has to vanish in a — 0 for all proper variations 

f ( V s c s | t=0 , c fl | fl=(U= o)ds +J J (&(v,V s q)c s , c n | a=0 ) + (V s V s q, V a q\ a=0 )dtds = 0. 

Unfortunately, we cannot yield any conclusions at this point, because c a ( 0, s, t) and 
V„r/(0, s, t ) cannot be chosen independently, since q is not any vector field along c 
but its image via the Square Root Velocity Function. Using the notations 

B[s, t) = &{q,V s q)c s {s,f), 

1 1 (V s V s o, c f ) 

D{s,t) = — =V s V s q(s,f) - -- ' c t {s,t), 

Vied 2 | c t \ bu 



and the lighter notation u{t\) t,,t2 = pj}’ t2 (u(t\)) to denote the parallel transport along 
a curve c of a vector field u, the following simple manipulations provide us with a 
solution 


<V s c s (s, 0), c fl (0,s,0)}ds + 


ff' 


(B{s, t),c a (0,s, r)> + (D{s, t),V t c a i 0, s,t))dtds 


f 

= [ (V s c s (s,0),c a (0,s,0))ds+ f f /B(s, t), c a (0, s,0) 0,t + f V f c (! (0,s,T) T ’ r dr\dfds 
Jo Jo Jo \ Jo / 

+ [ [ (D(s,t),v t c a (0,s,t))dtds 
Jo Jo 

= [ (v s c s (s,0)+ f B(s,t) t,0 dt,Ca(0,s,0)\ds+f f f (B(s, t) t,T ,V t c a (0, s,t)) 

Jo \ Jo / Jo Jo Jo 


drdfds 


(D(s, t),v r c a (0,s, t)> dfds 


+// 

= J o ( V s c s (s,0) + J o B(s,t) t ’°dt,Ca(0,s,0) 

+ Jo Jo ( D ^ S ’ t ' >+ J, B ( S ’ T ) T ' tdT ’ V tCa(0,S,t)^ 


dtds. 


Since the variations c a (0, s, 0) and V t c a (0, s, •) can be chosen independently for all s, 
we obtain the following characterization of the geodesics s >-► c(s, •) of ..M 


V s c s (s ; 

D(s ; 


, 0 ) + f 

Jo 


S&{q,V s q)c s {s, t) t,0 dt - 0, Vs 


V s g)c s (s,T) T,( dT = 0, Vt,s. 


Finally, using the definition of D(s, t), the geodesic equations are 
V s c s (s,0) = r(s, 0), Vs 
V s V s q(s, f) = || < 7 (s, f)|| (r(s,f) +rfs.f/), Vf,s 

where q is the SRV representation of c, the vector field r is given by 


(9a) 

(9b) 


r (s,t) = - J" &(q,V s q)c s (s,T) T ' r dT, 


and - (r, v) v with v = y 7~\ c > >' s the tangential component of r. 

6 Exponential map 

In this section, we describe an algorithm which allows us to compute the geodesic 
s >-► c(s, •) starting from a point c £ .M at speed u e T C M. This amounts to finding the 
optimal deformation of the curve c in the direction of the vector field u according to 
our metric. We initialize this path s >->• c(s, •) by setting c(0, ■) - c and c s (s,0) = u, and 



we propagate it using iterations of fixed step e > 0. The aim is, given c(s, ■) and c s (s, •), 
to deduce c(s + e, •) and c s (s + £,■)■ The first is obtained by following the exponential 
map on the manifold M 

c(5 + £,-)= exp c(v) ec s (s, ■), 

and the second requires the computation of the variation V 5 c 5 (s, •) 
c s (s + £,•)= [c 5 (s, ■) + eV s c s (s, •)] s,s+e , 


where once again, we use the notation w(s) s,s+e = P^’ s+e ( w(s )) for the parallel trans¬ 
port of a vector field 5 >-► w(s) along a curve 5 >-► c(s) in M. If we suppose that at time 
swe have c[s, •) and c s [s, •), then we also know Cf(s, •) and q(s,-) = -^==Cf(s, ■), as well 
as V f c s (s, •) and 


V s £/(s, ■) 


y s c t 

y/\ct\ 


(v) 


1 (VsCf.Cf) 

2 |cd 5/2 


C t (v)> 


( 10 ) 


using the fact that V 5 c f = V t c s . The variation V s c s (s, ■) can then be computed in the 
following way 


V s c s (s, t) = V s c s (5,0)°’ f + f [V s V 5 Cf(S,T)-5g(Cf,C s )C s (S,T)] T ' f dT 
Jo 

for all f e [0,1], where V s c 5 (s, 0) is given by equation (9a), the second order variation 
V s V s c t (s, •) is given by 


, „ (VfCc, c t ) 

V s V s c f = [Cf| 1 / 2 V S V S <7 + f/ V t c 5 

led 2 

3{V t c s ,c t ) 2 
2' 


(y stfi c t ) 


led 


3/2 


led 4 


|V/cd 2 

led 2 


c t , 


( 11 ) 


and V s V s q can be computed via equation (9b). 

Now, if we no longer have a continuous curve but a series of discrete observa¬ 
tions po,pi,...,p n £ M made at discrete times 0 = to < t\ < ... < t n = 1, then the 
optimal deformation of these points in the direction of a series of tangent vectors 
Mo, Mi,..., u„, where u k £ T Pk M for k = 0,..., n, can be computed by the following 
steps. 


Algorithm 1 (Discrete Exponential Map) 

Initialization : Set c(0, t k ) = p k and c s (0, t k ) = u k for Ic n. 

Heredity: Ifc[s, t k ) and c s ( 5 , t k ) are known for all k, then 


1. for all 0 < k < n - 1, set 


c t (s,t k ) 


V f c s {s,t k ) 


t 1 , lo gcfetjt) c($, ffc+t), 
l k +1 Ik 

-——- (Ci(s, t k+l ) tk+1 ’ tk - c s (s, t k )) , 
t k +1 - t k 


where log denotes the inverse of the exponential map on M, and compute q[s, t k ) = 
■^j=Ct(s, ffc) andV s q(s,t k ) using equation (To). 












2. Compute r(s, t k ) = E” =jt 2%{q, S/ s q)c s [s, t() te,tk for all 0 < k < n and set 

V s c s (s, to) = r(s,t 0 ), 

V s V s q(s,t k ) = ||< 7 (s, ffc)||(r(s, t k ) + r(s, ft) 11 ). 

3. Fork - 0,...,n—l, 

(a) compute^/ s SZ S c t {s, t k ) using equation UD- 

(b) computes/ s c s {s,tk + i) using 

S/ s c s [s,t k+ 1 ) = [ V s c s (s, t k ) + (ffc+i- t k ) (V s V s c f (s, t k ) - 3%(c s ,c t )c s {s,t k )) ] tk ’ tk+1 

4. Finally, for all 0 < k < «, set 


c(s + c, f fc ) = exp c(ijfjt) (ec s (s, t k )) 
c s ls + e, t k ) = [ c s (s, f fc ) +eV 5 c s (s, t k ) ] tk - tk + k 


7 Conclusion 

In the same way that the first-order Sobolev metric {2} on the space of plane curves 
can be obtained as the pullback of the L 2 -metric by the square root velocity func¬ 
tion (ED), our metric G can be obtained as the pullback of a natural metric G on 
the tangent bundle T Jl by the same SRVF. As such it is reparameterization invari¬ 
ant, and induces a Riemannian metric G on the shape space SP for which the fiber 
bundle projection is formally a Riemannian submersion. On the other hand, the spe¬ 
cial role that G gives to the starting points of the curves induces another formal fiber 
bundle structure, this time over the manifold M seen as the set of starting points of 
the curves, for which the projection is formally also a Riemannian submersion. The 
geodesic distance induced by G takes into account the distance between the origins 
of the curve and the L 2 -distance between the SRV representations, but instead of 
transporting the computations in a unique tangent plane as in (S) and I H i, we stay 
in the manifold. This should allow us to take into account a greater amount of infor¬ 
mation on its geometry. Finally, explicit equations can be obtained for the geodesics 
on the space Ji of parameterized curves for our metric G, and the exponential map 
can be iteratively computed. Future work will include some numerical simulations 
to illustrate the work of this paper. 
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